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Abstract 

We study a problem of optimal investment / consumption over an infinite horizon in 
a market consisting of two possibly correlated assets: one liquid and one illiquid. The 
liquid asset is observed and can be traded continuously, while the illiquid one can be 
traded only at discrete random times corresponding to the jumps of a Poisson process 
with intensity A, is observed at the trading dates, and is partially observed between 
two different trading dates. The problem is a nonstandard mixed discrete/continuous 
optimal control problem which we face by the dynamic programming approach. When 
the utility has a general form we prove that the value function is the unique viscosity 
solution of the HJB equation and, assuming sufficient regularity of the value function, 
we give a verification theorem that describes the optimal investment strategies for the 
illiquid asset. In the case of power utility, we prove the regularity of the value function 
needed to apply the verification theorem, providing the complete theoretical solution of 
the problem. This allows us to perform numerical simulation, so to analyze the impact 
of time illiquidity in this mixed market and how this impact is affected by the degree 
of observation. 
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1 Introduction 

Following the seminal works of Merton on portfolio management, a classical assumption in 
mathematical finance is to suppose that assets may be continuously traded by the agents 
operating in the market. However, this assumption is unrealistic in practice, especially in 
the case of less liquid markets, where investors cannot buy and sell assets immediately, and 
have to wait some time before being able to unwind a position. 

In the recent years, several works have studied the impact of this type of illiquidity 
on the investors. Rogers and Zane [23], Matsumoto [20] . Pham and Tankov [21] (see also 
[5[ [22]) consider an investment model where the discrete trading times are given by the 
jump times of a Poisson process with constant intensity A > 0. Bayraktar and Ludkovski 
[3] study a portfolio liquidation problem in a similar context. 

The aforementioned works focus on an agent investing exclusively in an illiquid asset. 
However, in practice it is common to have several correlated tradable assets with different 
liquidity. For instance an index fund over some given financial market will be usually 
much more liquid than the individual tracked assets, while sharing a positive correlation 
with those assets. An investor in this market will then have the possibility of hedging his 
exposure in the less liquid assets by investing in the index and rebalancing his position 
frequently. 
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To our knowledge only few papers consider the case of a market composed by two 
(possibly correlated) assets, one liquid and one illiquid. This is the case of Longstaff |19j . 
who analyzes a two agents portfolio problems in a market composed by a liquid asset 
and another asset that becomes non tradable only for a given time period. Moreover, 
Schwartz and Tebaldi [M] consider a market composed by a liquid asset that can be traded 
continuously, and an illiquid asset that cannot be traded and is liquidated at a terminal 
date. Finally, we have to mention the recent paper by Ang, Papanikolaou and Westerfield pQ 
where the authors, in the standard framework of maximizing an infinite horizon discounted 
power utility from consumption, take a less restrictive point of view on the tradability of 
the illiquid asset assuming, as in [5j [121 ISHl EH [23], that the illiquid asset may be traded 
at discrete random times. 

Following the line of these papers, here we also consider a market composed by a 
liquid asset and an illiquid one. In particular, as in pQ the illiquid asset can be traded 
at some discrete random dates. The main novelties of our paper in the context of mixed 
liquid/illiquid market case are the following: on one hand, we consider the case of incomplete 
observation on the illiquid asset price between trading dates; on the other hand, we provide 
a complete theoretical framework that allows to cover also more general form of utility 
functions, not restricted to CRRA type. 

To be more precise, we study a problem of optimal investment/consumption over an 
infinite horizon in a market consisting of a liquid and an illiquid asset. The liquid asset is 
continuously observed and can be continuously traded. The illiquid asset is correlated with 
the liquid one with correlation parameter p E (—1,1) and can be traded only at discrete 
random times corresponding to the jumps of a Poisson process with intensity A > 0. About 
the observation, we assume that the illiquid asset can be observed at the trading dates 
(as in [5j [I2j [20l [21]), but we introduce a new feature in the model - with respect to the 
aforementioned literature - allowing to consider the possibility of a partial information 
between trading dates. We introduce a parameter 7 £ [0, 1] measuring the observation of 
the illiquid asset between two trading dates. The limit cases for this parameter, i.e. 7 = 
and 7 = 1, correspond respectively to the no observation case (as it is done in [51 [2T| 122] . 
but only for the single asset case) and to the full observation case (as in [TJ [23| 124]). In 
this sense our model is more general and flexible than the other ones proposed by the 
literature: first, it allows to put together the presence of two correlated liquid/illiquid 
assets and the incomplete observation of the illiquid one; second, to merge two extreme 
(full vs. none observation) cases and consider intermediate situations where the agent has 
a partial information about the state of the illiquid asset between two trading dates (see 
Remark E2]). 

The mathematical problem is, as expected, a nonstandard mixed discrete/continuous 
optimal control problem, which is in our case more difficult to tackle than in the afore- 
mentioned papers due to the presence of the features described above. By means of a 
suitable use of Dynamic Programming (DP) (extending what is done in |21j). we show that 
the stochastic control problem between trading times can be written as an infinite horizon 
stochastic time-inhomogeneous control problem. Then we apply the usual machinery of 
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DP for such kind of problems and, using some results of [7] I . we characterize the value 
function V of this auxiliary problem as the unique viscosity solution of a Hamilton-Jacobi- 
Bellman equation. At this stage, the viscosity characterization only allows to provide the 
optimal feedback allocation in the illiquid asset. In order to go further in the solution and 
characterize the optimal feedback allocation in the liquid asset as well as the optimal feed- 
back consumption strategy, we need to prove a regularity result for V. Such (nonstandard) 
regularity result is provided in the special case of power utilitjH and allows to give a full 
theoretical solution to the problem. This solution is made implementable by a numerical 
scheme and numerical results are then provided and discussed for different values of the 
relevant parameters 7, A, p. 

The plan of the paper is as follows. Section 2 describes the market model and formulates 
the investment/consumption problem. In Section 3 we show how, by a suitable dynamic 
programming principle, the problem can be reduced to a standard continuous time stochas- 
tic control problem; then we present some useful properties satisfied by the value functions 
- the original one and the auxiliary one - and prove an analytical characterization of them 
by means of viscosity solutions; finally, we provide the characterization of the optimal in- 
vestment in the illiquid asset. In Section 4 we solve completely the problem in the case of 
power utility and provide an iterative scheme to solve the problem numerically. Finally, 
Section 5 is devoted to the discussion of the numerical results we obtained. 



2 Model and optimization problem 

In this section we present the model and the optimization problem we deal with. Let us con- 
sider a complete filtered probability space (fi, T, (7i)(>o,P) satisfying the usual conditions, 
on which are defined: 

- A Poisson process (N t ) t >o, with intensity A > 0. We denote by (rk)k>l its jump 
times; moreover we set to = 0. 

- Two independent standard Brownian motions (Bt)t>o, (Wj^o, independent also of 
the Poisson process (N t )t>o- 

2.1 Market model 

The market model we consider on the probability space above consists of two risky assets 
with correlation p € (— 1, 1): 

- A liquid risky asset that can be traded continuously; it is described by a stochastic 
process denoted by L t whose dynamics is 

dL t = L t (b L dt + a L dW t ), (1) 

1 Actually, in [7] these results are proved for the case 7 = 0. However, their extension to the general case 
7 € [0, 1] is straightforward, see Subsection 13. II 

2 Our assumption on the utility function covers only the case of positive power, differently from [T|. 
However the methods used here can be employed by suitable modifications to cover the case of negative 
power as well (see also Remark 12. 4|l . 
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where &l, gl > 0. 

- An illiquid risky asset that can only be traded at the trading times r^; it is described 
by a stochastic process denoted by It, whose dynamics is 

dl t = h (bidt + en {pdW t + P 2 dB t fj , (2) 

where bi, ai > 0. 

Without loss of generality we assume Lq = Iq = 1. We also suppose that in the market 
is present a riskless asset with deterministic dynamics. Just for simplicity we assume that 
the interest rate of this asset is constant and equal to 0. 

Remark 2.1 If the riskless asset interest rate is not 0, one needs just to add an extra term 
in all the equations. Moreover, in the special case of power utility treated in Section^ the 
assumption that the riskless rate return is can be done without loss of generality, as in 
this case the interest rate can be discarded in the discount factor of the objective functional 
(the constant /3 m (jSJ) below) by a suitable change of variables (see Remark 2, p. 189, in 

w- 

2.2 Information 

The information setting we consider is the following. We assume that: 

- the liquid asset L is continuously observed; 

- the illiquid asset / is observed at the trading random times (Tk)k>o', 

- the illiquid asset / is only partially observed in the time interval (rfc,Tfc+i). 

In order to deal with and make formal the last issue, we suppose that the Brownian motion 
Bf can be split as 

B t = 7 S t (1) + s/T=i*dBi 2 \ 7 G [0, 1], 

where B^ l \B^ are mutually independent Brownian motions also independent of W,N, 
and we assume that B^ is observed and B^ is unobserved. 

Let (Mt)t>o, (Wt)t>o, {b\ )t>o be the nitrations generated respectively by N, W, B^. 
Define the u-algebra 

It 

Moreover define the filtration 

G° := {Gt)t>o; Qt 
The observation filtration we consider is 

G = (Gt)t>o; Qt = Qt V ff(P-null sets). 
This means that at time t we have: 



a 



{r fe <t}> 



k > 0) , t>0. 



Af t V l t V W t V 6 f (1) = a(r k , I Tk ; r k < t) V W t V . 
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- a full information on the past of the liquid asset up to time t; 

- a full information on the trading dates of the illiquid assets occurred before t and on 
the values of the illiquid asset at such trading dates; 

- a partial information (in the sense described above) on the value of the illiquid asset 
at time t. 

The parameter 7 measures how much information on / is available in the random intervals 
(rfc,Tfc_|_i). The limit cases are: 

-7 = 0, which corresponds to have no information on the illiquid asset in the interval 
(Tfe,Tfe+i) after the last information on its value at time r^; this case fits the model 
described in [22]; 

-7 = 1, which corresponds to the full information case, recovering the setting of pQ. 

Remark 2.2 Roughly speaking, the main idea behind our model is that the price of the 
illiquid asset is observed exactly at the trading times (rk)k>o (as in \2Q \2^), while at 
different times t G (Tfc,Tfc +1 ) the agent can observe a process I^ 1 satisfying the equation 

dl[ l) = ll 1) {b I dt + a I { P dW t + ^/l^ 1 dB^ ) )), jW = I Tk , 

and 

r r(l) r(2) 

l t - H - H > 

where 

is an unobserved noise on our knowledge of I. The idea behind this is that, while the price 
of I is clearly observed when it is bought/sold, between two trading dates we have only a 
partial knowledge of it. This partial knowledge is represented by the process and differs 
by the "real price " I by a factor 1^ . 

2.3 Trading/consumption strategies and wealth dynamics 

In the setting above, we define the set of admissible trading/consumption strategies in the 
following way. Consider all the triplets of processes (c, tt, a) such that 

(hi) c = (ct)t>o is a continuous-time nonnegative process (^t)t>o-predictable and with 
locally integrable trajectories; q represents the consumption rate at time t; 

(h2) 7r = (7Tt)(>o is a continuous-time process (C?t)t>o-predictable and with locally square 
integrable trajectories; nt represents the amount of money invested in the liquid asset 
at time t; 

(h3) a = (afc)fceNj is a discrete process where is ^-measurable; represents the 
amount of money invested in the illiquid asset in the interval (r^, T/t-j-i]. 
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Given an initial wealth r > and a triplet (q, nt, ot k ) satisfying the requirements (hl)-(h3) 
above, we can consider the process Rt representing the wealth associated to such strategy. 
Its dynamics can be defined by recursion on k > by 

R = r, (3) 
R t = R Tk + J {-K s {b L ds + a L dW s ) - c s ds) + a k (^j- - l^j , t G (r fc , r k+1 }. (4) 

We observe that the process R is in general not G-adapted (unless 7 = 1), as / is not. 

We can define also a liquid part Xt (observable, therefore G-adapted]f| and an illiquid 
one At (partially observable) of the wealth Rt- They are defined in the intervals [r k , T k+ i) 

X t = R Tk -a k + / (ir s (b L ds + a L dW s ) - c s ds), (5) 



A t = a k -±. (6) 

Obviously we have 

Rt = Xt + At, Vt>0. 

Observe that the process R is continuous, while the processes X, A are not, due to the 
rebalancing. We also observe that at time r k , for any k, the process R does not depend on 
the value of a k , while the process X, A do. 

As a class of admissible controls we consider the triplets of processes (c, ir, a) satisfying 
the measurability and integrability conditions above and such that the corresponding wealth 
process Rt is nonnegative (no-bankruptcy constraint). One can see without big difficulty 
that this requirement is equivalent to require that both the liquid and the illiquid wealth 
have to be nonnegative at each time, i.e. that Xt > 0, At > for every t > 0. So, the 
admissibility of a strategy (c, w, a) is equivalent to 

< a k < R Tk , Vk > 0, (7) 

t 

(ir s (b L ds + a L dW s ) - c s ds) < R Tk - a k , Vi G [r fc ,r fc+1 ), V/c > 0. 

k 

The class of admissible controls depends on the initial wealth Rq = r. We denote this class 
by A(r) noticing that it is not empty for every r > 0, as the null strategy (c, tt, a) = (0, 0, 0) 
belongs to it for every r > 0. 



2.4 Optimization problem 

Let Rq = r. The optimization problem consists in maximizing over the set of admissible 
strategies A(r) the expected discounted utility from consumption over an infinite horizon. 

3 Notice that X t is the sum of the money invested in the liquid asset and of the money held in the bank 
account (which has interest rate by our assumption) . 
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In other terms, chosen a utility function U and a discount factor (3 > 0, the optimization 
problem we consider is the mixed discrete/continuous stochastic control problem 



Maximize E 



e- ps U{c s )ds 



o 



over (c, 7r, a) G „4(r). (8) 



Assumption 2.3 T/ie preference of the agent are described by a utility function 

U : R+ — > M 

which is continuous, nondecreasing, concave, and such that U(0) = 0. Moreover we assume 
the following growth condition on U: there exist constants K{j > and p € (0, 1) such that 

U(c) < K v -. (9) 
p 

Remark 2.4 In the applications one is usually interested to work with power utility func- 
tions of the form 

c p 

U(c) = — , pe(-oo,l), 
p 

with the usual agreement that U(c) = logc when p = 0. However, Assumption \2.3\ includes 
only the case p £ (0,1). On one hand the case of negative exponent is interesting, as it 
seems to be even more realistic from the point of view of the agents' behavior (see f^); on 
the other hand, to be extensively treated it would need much more space. For this reason we 
will work with Assumption \2.3[ Nevertheless we stress that the case p < can be treated by 
the same techniques by straightforward modifications. This case is treated in J7}/ under full 
observation. 

We observe that, due to Assumption 12.31 on U, the Legendre tranform of U 

U(w) = sup {U(c) — cw}, w > 0, 

c>0 

is finite, nonincreasing and convex. Moreover, the growth condition © yields the following 
growth condition for U: there exists > such that 

U(w) < K v w~^. (10) 
According to Remark 12.61 b elow . we assume the following: 
Assumption 2.5 We assume that 

P > k p , (11) 

where 

f pO- ~ p) 1 

k p := sup {piuLbL + uib!) (u 2 L a 2 L + ujaj +2pu L uia L cr I )\ . (12) 

u L eR, Ul £[o,i] I 1 ) 
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Remark 2.6 The assumption on f3 is related to the investment/ consumption problem with 
the same assets but in a liquid market. Let p E (0, 1) and consider an agent with initial 
wealth r, consuming at rate Ct and investing in Lt and It continuously with respective 
proportions u\ and u\ and under the constraint that u\ G [0, 1] . Suppose, moreover, that 
the preferences of the agent are represented by the utility function U^ p \c) = c p /p, with 
p £ (0, 1). Let us denote by AMert{ r ) the set of strategies keeping the wealth nonnegative 
and define the value function 



(p) 



Alert 



(r) = sup E 

(u L ,U ! ,c)eA M ert(r) 



-^U {p \c t )dt 



(13) 



This is a constrained Merton problem which dominates our problem, in the sense that 
V$ Ar) is higher than the optimal value of our problem, up to the multiplicative constant 
Kjj of 0. One can see (for instance by solving the HJB equation) that V^ ert is finite if 
and only if (fTTj) is satisfied and that in this case 



Mert 



(r) 



1-P 



i-p 



(14) 



Therefore, condition (|lip guarantees together with §§§ finiteness for our problem too. 

Further note that the constrained liquid investment/consumption problem described above 
can always be reduced to the case where the two assets are independent, because 



dX t = X t (uf 



dL t , jdlt 



u t + 



o-l 



-iii 



dLt 
L t 



+ n 



■ dJ t 



where J is the process defined below in (I21j) ( taking 7 = 0), and the problem is equivalent 
to an agent investing in L and J, with the same constraint for the proportion invested in 
L. However, this reduction does not work for the illiquid problem that we consider: neither 
the observation constraint (the integrand in L being G-adapted) nor the trading constraint 
(the amount held in the illiquid asset being constant between and Tk+i) are preserved by 
this transformation. 



3 Dynamic Programming 

We denote the value function of the optimal stochastic control problem ([8|) by V: 



V(r) 



sup E 

(c,7T,a)eA(r) 



-P s U{c s )ds 



r > 0. 



(15) 



Proposition 3.1 V is everywhere finite, concave, p-Holder continuous and nondecreasing. 
Moreover 



V(r) < Kyr p , for some Ky > 0. 



(16) 



9 



Proof. As we have already observed in Remark 12 .61 finiteness and (|16p follow from ([9]) and 
(|lip . by comparing with a constrained Merton problem. 

Concavity of V comes from concavity of U and linearity of the state equation by stan- 
dards arguments. Also monotonicity is consequence of standard arguments due to mono- 
tonicity of U. Finally p-Holder continuity follows from concavity and monotonicity of V 
and from (|16p . □ 



Following [21], we state a suitable Dynamic Programming Principle (DPP) to reduce our 
mixed discrete/continuous problem to a standard one between two trading times. 



Proposition 3.2 (DPP) We have the following equality: 



V(r) 



sup E 



e-P s U{c s )ds + e- liT1 V (R T1 ) 



(17) 



Proof. The proof is long and technical, but similar to the one in [22] and we omit it for 
brevity. However, we note that, unlike in [22], there is some additional random information 
between and n brought by W, so that the "shifting" procedure is slightly more technical 
to achieve. We refer for instance to the proof of the DPP provided in Appendix B of [12 , 
where this shifting procedure is performed in a similar framework. □ 



Now we will use this DPP to relate our original problem to a standard continuous-time 
control problem. For each x > 0, let Ao(x) be the set of couples of stochastic processes 
(c s ,7r s ) s >o such that 

- (c s ) s >o is (W a V £>i^) s >o-predictable, nonnegative and has locally integrable trajec- 
tories; 

- (vr s ) s >o is (Ws V 2?s ) s >o-predictable and has locally square-integrable trajectories; 

- x + Jq \—c s ds + iT s (bLds + (JLdWs)) > 0, for every T > 0. 
By Lemma |A.1| we obtain that ()17p may actually be rewritten as 

e-P'U(c a )ds + e- pn V (R^ 



V(r) 



sup sup E 

0<a<r (c,ir)<=Ao{r-a) 



(18) 



We want to rewrite in a suitable way the inner optimization problem in (|18p . i.e. 



sup E 

(c,Tr)£Ao(r— a) 



f 1 e-P s U(c s )ds + e~P T1 V (R 7 
Jo 



(19) 



First define (see Remark 13. 3f i) for explanations on this choice) 



by = 7 6/ + (1 - 7 ) . 

O-L 



h = (l-7 2 )(^-^), 



(20) 
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and given x,y > 0, (c, 7r) E ^4o(^)) define the processes J, X x ' c,7r , Y y as solutions to the 
SDEs 



dJt = Jt (bjdt + a iy /l - pVl - l 2 dBf ') . ^0 = 1, 

cL¥ s = -C s ds + TT s Q)Lds + (TLdW s ), X = x, 

dY s = Y s (b Y dt + <r 7 (pdW s + Vl " P 2 7^ (1) )) , 



o = y- 



(21) 
(22) 
(23) 



Then we have, for every t £ [0, t\) 



X t = X. 



A t = Y t a ° ■ J t . 



Since t\ is independent of (Wt)t>o, (B^)t>o, (B^)t>o and has distribution S (A), while c, 
J, X*&«,YV are (Woo V V B^-measurable, we have 



E 



e-P s U(c s )ds + e-^ViRr,) 



/3ni 



Xe 



e-P s U(c s )ds + e- /3 V(A > [- a ' c '" + J t ■ Y t a )j dt 

f'OO f'OO f'OO 

/ e-P'Uic) / Xe- Xt dtds+ / Ae-( A+ ^V(Xr a ' C,7r + Jt ■ Y t a )dt 

JO Js Jo 



U(c t ) + \V{X r - a ^ + J t ■ Y t a ) dt 



where in the second equality we have used Fubini's theorem. 

On the other hand, since J is independent of (Wt)t>o and (B^ )t>o, while c, X x ' c,7T ,Y y 
are (Woo V £>co^)-predictable, conditioning the equality above with respect to Woo V Bdo we 
get 



E 
E 



e- ps U{c s )ds + e-P Tl V(R 



Tl - 



-(/3+A)t 



U(c t ) + \V(X r t - a ^ + J t ■ Y t a ) )dt 



WocVfiW 



(c/( Ct ) + AG[y](t,X 



r— a,c,7r v n 



^ a dt 



where 

G[Vl(t,x,y) := E[^(x + yJ t )]. 
In conclusion, we may rewrite (|19|) as 



sup E 

(c,7r)e»4o( ? '— a ) 



-(/3+A)t 



c/( Ct ) + AG[y](t,xr^,y t a )Ut 



(24) 



(25) 



The new form of the problem involves (|24p . It is useful to define G as a linear operator 
from the space Ai i (R + ; R) of measurable functions with at most linear growth to the space 
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of measurable functions M. (M+ ; K) and stress its dependence on 7 (it depends on 7 through 
J). So we define 

G 7 : 7Wi(M + ;M) — ► M(K;j.;R) (26) 
i> .— > G 7 [V#,*,y) := E + yJ t )] . (27) 

Useful properties of G 7 are listed in Proposition IA.2I in Appendix. 

The problem (|25p is a continuous non autonomous stochastic control problem over 
an infinite horizon that we call auxiliary problem. One can apply a standard dynamic 
programming approach to this problem. To do that, as usual we define the same problem 
for generic initial data t, x, y. For each t, x > 0, let At(x) be the set of couples of stochastic 
processes (c s ,7r s ) s > t such that 

- (c s ) s >t is (Ws V ^ 1 ^) s >t-predictable, nonnegative and has locally integrable trajecto- 
ries; 

- (^s)s>t is (W s V £>i^),s>i-predictable and has locally square-integrable trajectories; 

- x + ff(- Cs ds + TT s (b L ds + a L dW s )) > 0, VT > t. 



Let x, y > 0. Given (c, n) € At(x), let (X s ' x ' c ' n ) s >o, (Y s ' y ) s >o be the solutions to the SDE 

dX s = -c s ds + ir s (b L ds + a L dW s ), X t = x, (28) 
dY s = Y s (b Y dt + a^pdWs + Vl " P 2 7^ (1) )) , Y t = y. (29) 



Define the value function 
V(t, x, y) = sup E 

(c,ir)&At(x) 



oo 

e -(H+X){s-t) 



(U(c s ) + XG^V^s,^^,?^ ds 



• (30) 



Associating to every locally bounded function v on the function Hv defined on M + by 

[Hv](r) = sup v(0, r — a, a), 

0<a<r 

by the arguments above we may rewrite the original value function (|15|) as (see also (|18|) 
and (J30J)) 

V(r) = [HV](r). (31) 



The problems (|30|) - ([3T|) are coupled since V is defined in terms of V in (|30|) . and, viceversa, 
V is expressed in terms of V by (|31l) . 



Remark 3.3 ^ H^e explain our particular choice for the drifts by and bj in (|20p . -Firsi 
of all we clarify that what we need for our argument is a couple of processes (Y, J) such 
that: (i) Y ■ J = A on [0, n), where A is defined in ©; (ii) Y is G-adapted; (Hi) J is 
independent ofG. 



4 Note that while we take the state variable X, Y starting at t, the process J remains the same as in ()21[) 
and enters in G 1 making the auxiliary problem nonautonomous. 
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Therefore, it is natural to consider the processes driven by SDE's (|2ip and (|23j) . with 
by, bj that can be chosen freely under the constraint bj + by = &/. Define the constants 

kL,Y,p ■= sup \ p(u L b L + Uyby) 

u L £R,u Y £[0,l] 1 



p(l-p) 



[u\a\ + u 2 Y aj(p 2 + 7 2 (1 - p 2 )) + 2pu L u Y o- L a I )^, (32) 



fc JlP := sup Sp bjUj -^—Pla 2 (l-p 2 )(l- 7 2 ) U 2 ). (33) 

These constants naturally appear respectively in Lemma \A.3\ and Proposition \A.2V (v). Com- 
bining these two results with (|16p . one gets an estimate on the growth ofV (precisely esti- 
mate (|37p below) under the condition that f3 > fc^yjj + fcjjj. in Lemma \A.J\ in the Appendix 
it is proved that, for our choice ofby,bj, we have 

kL,Y,p + k J,p = k p , (34) 

which is the minimum possible value of ki,Y,p + kj p . So, our choice of drifts allows then to 
treat the auxiliary problem without further restriction on (5 other than A s sumption 1 2. 51 

(%%) We notice that the auxiliary problem (|3Q|) is not autonomous due to the dependence 
of G 7 [y] on time. However, in the case of full observation (j = 1), one has ,7=1, and 
G l \V](t,x,y) = V(x+y). In this case, consistently with Jl]/ ; we get an autonomous problem^ 



sup E 



e -C0+A). (V(c s ) + XV (x^'C + Y s °' y ^ ds 



Therefore, time inhomogeneity of our auxiliary problem is due to the lack of full information. 
It is also worth to outline here that in (|3Up we take the discount e~^ +A ^ s_ *) in place of the 
more natural e~^ + ^ s , as it allows to get rid of the exponential terms in the HJB equation. 

(Hi) The value function V is the analogue of the value function of JH \2^j. Indeed, in 
Hi\2^, the initial time is not supposed to be a trading time (of course this can be done also by 
us without loss of generality), so the initial endowments x,y in the liquid and illiquid asset 
cannot be rebalanced optimally at t = - in other terms, the analogue of the optimization 
(13 1|) is not performed at t = in fJl \2J$ . 

3.1 HJB equation and viscosity characterization of V 

In this section we characterize V as unique viscosity solution of an associated Hamilton- 
Jacobi-Bellman (HJB) equation. We start stating some qualitative properties of V. The 
proof can be found in [7] in the case 7 = 0. In the general case the proof proceeds exactly 
in the same way and we omit it for brevity. 



5 With respect to pQ, our dynamic programming approach to the problem is different. Indeed, apart 
the fact that we deal with general utility functions (which is not relevant for that), our approach seems to 
be the only possible to deal with the issue of possible partial information (7 < 1). Hence, the differential 
problem we get is different from the one derived in [T], also when our control problem coincides with the 
control problem of [1] (7 = 1 and power utility). Nevertheless, as suggested by the intuition, our differential 
problem must be autonomous as well as the one of [T], when the two control problems coincide. 
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Proposition 3.4 V(t, ■) is concave with repect to (x,y) and nondecreasing with respect to 
x and y for every t > 0. Moreover it satisfies the boundary condition 



v{t,o,y) = e 



oo 



Vt > 0, Vy > 0. (35) 



In particular, since by Assumption \2.3\ it is U(0) = 0, due to Proposition \A.^( v) and (fl~6j) . 
roe have 

V(t,0,0) = 0, Vt>0. (36) 
Finally, V is continuous on and satisfies, for some > 0, the growth condition 

< V(t,x,y) < K 9 e k ^\x + yY, y(t,x,y) G . (37) 



By standard arguments of stochastic control (see e.g. |25} Ch. 4]), we can associate to 
V an HJB equation, which in this case reads as^| 

- v t + (P + X)v - \G[Uv] - sup H cv (y,Dr x }V )V,D? X .v; c,ir) = 0, (38) 

c>0,7reK 

where for (y, q, Q) 6 R + x K 2 x £2 (where 1S2 denotes the space of symmetric 2x2 matrices), 
c > 0, 7T € M, the function H cv is defined by 

H C v(y,q,Q]c,7i:) 

2 2 2 

= 17(c) + (irb L - c)qi + b Y yq 2 + ^~Qn + *7WL2/Qi2 + {p 2 + 7 2 (1 - p 2 ))^y 2 Q22- 
Note that sup c>0 t€R H cv (y, q, Q; c, ir) is finite if q± > 0, Qn < 0, in which case we have 



sup H cv (y,q,Q;c,ir) 

c>0,7reJR 



(b L qi+ pCTLO-iyQu) , / 2 , 2/i 2\N°7 2/ 

tttt^ + b yyi2 + {p + 7 (i - p ))-±y < 



We are going to characterize V as unique constrained viscosity solution to ([38]) according 
to the following definition. 

Definition 3.5 (1) An upper-semicontinuous (resp. lower-semicontinuous) function v is a 
viscosity subsolution (resp. supersolution) to ([38]) at (to,xo,yo) £ if 

-cp t (t ,xo,yo) + W + X)v(t ,x ,yo) - XG[Hv](t ,x ,yo) 
- sup H cv (y ,D( xy )ip(to,x ,yo),D? x ) ip(to,xo,yoy, c,7r) < (resp. > 0) 

c>0,7reK 



6 Note that the HJB equation associated to the auxiliary problem would have as third term — AG[V] and 
then should be coupled with (|3ip . This is the approached followed by [2 1 1 1221 . Here we follow a different 
approach by inserting directly (|31[1 , which yields a nonlocal term in the equation. 
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for any function ip £ C 1,2 (IR^_;M) such that y?(to, x$, yo) = v(to, xo,yo) an d <p > v (resp. <) 
on [to, to + e) x 0, for some neighboorhood O C of (xo, yo) and e > 0. 

(2) We say that a continuous function v is a constrained viscosity solution to (|38p if it 
is a subsolution on a supersolution on M + x (0, +oo) x M + and satisfies the boundary 
condition 

r r°° i 
v(t,0,y) = E^jf e-( /3+A ^ s -^AG 7 [^5](s,0,y^)dsJ, Vt > 0, Vy > 0. (39) 

Remark 3.6 (i) To simplify some proofs later on (in particular for regularity properties 
of the solution), we have chosen a definition of viscosity solution that slightly differs from 
the usual one. Indeed, in the standard definition of viscosity solution for parabolic PDEs 
(e.g. W), the test functions ip should be above v in a neighborhood of (to, xo, Do), while in 
our definition we are only concerned with their position for t > to- We notice that our 
definition is in principle more restrictive than the usual one, as we are enlarging the set of 
test functions. However, as soon as comparison is satisfied for the standard definition the 
two definitions are actually equivalent, see 113)1 for results in this direction. 

(ii) The concept of constrained viscosity solution we use naturally comes from the 
stochastic control problem. The boundaries {x = 0, y > 0} and {x > 0, y = 0} are 
both absorbing for the control problem (in the sense that starting from these boundaries, 
the trajectories of the control problem remain therein), but they have different features. 
Indeed starting from the boundary {x > 0, y = 0} the control problem degenerates in a 
one dimensional control problem; the associated HJB equation is nothing else but our HJB 
equation restricted to this boundary and this is why we require viscosity sub- and superso- 
lution properties to the value function at this boundary. Instead starting from the boundary 
{x = 0, y > 0} there is no control problem (since At(0) = {(0, 0)}) and the natural condi- 
tion to impose is a Dirichlet boundary condition. 

Theorem 3.7 V is the unique constrained viscosity solution to (|38p satisfying the growth 
condition (|37p . 

Proof. The proof that V is a viscosity subsolution on and a viscosity supersolution on 
M + x (0, +oo) 2 follows exactly the same arguments of, e.g., [251 Ch. 4], once we note that 
such arguments look at test functions only for t > to- The Dirichlet boundary condition 
(|39l) is verified due to (|35|) and (|3T|) . The growth condition ([37]) has been already proved. It 
remains to show that V is a supersolution when y = 0. In this case, as noticed in Remark 
I3.6( ii). the control problem degenerates in a one-dimensional problem and again standard 
arguments apply to this control problem, giving the viscosity supersolution property. 

Uniqueness is consequence of the comparison principle Proposition 13.81 below, whose 
proof can be found in the case of standard definition of viscosity solution in [7] when 7 = 0. 
In the general case there is no special change to perform in the proof, so we omit it. □ 



Proposition 3.8 Let v\ (resp. V2) be a viscosity subsolution (resp. supersolution) to (|38p 
on x (0, +00) x R_|_. Assume that v\, V2 satisfy the growth condition (|37p . and the 
boundary condition 

r f°° 1 

vi(t,0,y) < E / e-^ +A ^ s- *>AG[?£t;i](s,0,Y^ I ')ds (40) 
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(resp. > for v-i)- Then v\ < V2 on R^_. 
3.2 Optimal policies 

In this section we show how to use the results obtained in the previous sections to construct 
optimal policies for the original problem. Unfortunately, since we do not know a priori if 
the value function V is smooth, we cannot write the optimal feedback for the controls (c, ir). 
However V can be computed numerically as viscosity solution of (|38p (see Section T4.6P and 
then an optimal allocation (a* k )k>o for the illiquid asset can be derived. Due to (|3T|) . at 
time t = an optimal allocation in the illiquid asset must be such that 

«o £ argmax < o < r V(0, r — a, a). (41) 

We notice that we cannot exclude that «q = at this stage. 

The expression of the optimal allocation in the illiquid asset (]41[) can be generalized to 
the random trading dates Tfe, for k > 1. At such trading dates, an optimal allocation must 
be such that 

a\ G argmaxo^^y^,^ - a, a). (42) 

This fact comes from the Markov property of our controlled system, whose proof we omit 
for brevity. 

Assuming sufficient regularitj0 for V we can prove a verification type result providing the 
optimal feedback for the controls (c,tt). 

Hereafter, by C l,k we will denote the class of functions which are once differentiable with 
respect to the time variable and /c-times differentiable with respect to the space variables, 
with continuous derivatives. 



Theorem 3.9 Suppose that 

Ui) V € C 1,2 (IR + x (0,+oo) 2 ;M), 

V"( v ,0) G C^ 2 (R + x (0,+oo);M). 

Suppose that there exist measurable feedback maps C* , n* such that 



(43) 



■ e argmax c > \U(c) - cV x (t,x,y)j , if x > 0, 
u (t, x, y) < 

:0, ifx = 0, 



U*(t,x,y) 

G argmax wm lirb L V x (t,x,y) + ^-V xx (t,x,y) + Kpoio L yV xy (t, x, y)j , if x > 0, 
= 0, ?/x = 0, 

(45) 

7 Such regularity will be proved in the next section in the case of power utility. 
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and that the closed loop equation 



dX s 
X t 



C*(s,X.,Y;> v )ds + U*(s,X s ,Y s l ^){b L ds + a L dW s ), 



(46) 



admits a unique nonnegative solution X*' t,x . Define the feedback strategies 
c* s = C*(s,X*/' x ,Y s ^), n* s = u*(s,x;^,y s ^). 

Then (c*,7r*) G At(x) and is an optimal control for the auxiliary problem (|3U|) , 



(47) 



Proof. First of all we prove that (c*, ir*) £ A t (x). Note that the fact that the SDE (gBJ has 
a well-defined solution implies that c* and 7r* satisfy the required integrability conditions, 
and moreover, by uniqueness of solutions of (|46p 



X 



t,X,C ,7T 



(48) 



concluding the proof of admissibility. 

To prove optimality we distinguish two cases: y > and y = 0. 

Case y > 0. In this case we have Ys' y > a.s. for every s > t. Therefore, due to the 
assumption of regularity (f4"3"|) -(i). we may apply Dynkin's formula to 



-(/3+A)0-t) 



V(s,X*/' x ,Y^) 



in the interval [t, r AT) for all T > t, where r = inf{s > t \ Xt' t,x = 0}. Since V solves (|38p 
in R + x (0, +oo) 2 in classical sense, by definition of C*,H* we get, arguing as in standard 
verification theorems, 



V(t,x,y)-E 



»tAT 



E 



-(/3+\)(s-t) 



UicD + XG^V^X*^,?^)) ds 



Splitting on the sets At = {t < T} and A T , we write 



V(t,x,y)-E 



1a t e 



-(/3+A)r 



E 



1a 7 



V(t, 0, Y T ^) + l^c e~^ T V(T, X*/> x ,Y*> y ) 
(P+\)(s-t) fjj^ + AGT[y]( s ,x s * Ax ,y^)) ds 

T 



+ 1a% 



-(/3+\)(s-t) 



U(c* s ) + XG^[V}(s,X*/' x ,Y^)) ds 



(49) 



By (|35p and using the Markov property of Y and the fact that At is Q T measurable, 



i AT v(T,o,Yt>y) 



E 



E 



1a 2 



Qr 



'^ +x ^ s - T) XG' r [V}(s,0,Y s t ' y )ds 



Qr 



(50) 
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Noting that X*' t,x = from r on and moving the term corresponding to (|50p to the right 
hand side in (09]), we get 



V(t,x,y)-E 



1a° e" 



E 



E 



-(p+\)(s-t) 



U(c* a ) + X(P[V](s,X*/' x ,Yly))ds 



• (51) 



Now we take T — > oo in (|5ip . The second term of the left hand side converges to 
by dominated convergence due to (fTTJ) . (fMj) . (|37|) and (J98J) . On the other hand, since 
J{t, x; c* , ir*) < oo, we have almost surely 



Hence, 



s l [T , +oo) ( S )e-^ +A )^^( C :) + AG^y]( S ,X^,y s ^) 

is almost surely dominated by a function which is integrable on [t, +oo) and goes to when 
T — > oo. This shows, by dominated convergence that almost surely 



-Q3+A)(s-t) 



[/•(<$) + AGT[F](s,x;^,y^) ) ^ 



0. 



Then by dominated convergence on Q we get that the last term of the right hand side of 
([5TD goes to 0. 

Case y = 0. In this case iV 2 ' = a.s. for every s > 0. This means that the boundary 
{x > 0, y = 0} is a trap for the control problem, which indeed lives in a one-dimensional 
space. The value function of this control problem is V(-, - ,0) and the HJB equation is still 
(|38p restricted on the half-line {x > 0, y = 0} (all the terms containing the derivatives with 
respect to y indeed disappear). Due to the regularity assumption (|35])-(ii). a verification 
argument similar (indeed easier) to the one above yields the claim. □ 

In the next section we will be able to show that the assumptions of the theorem above 
are fulfilled in the case of power utility. We can put together (02]) and Theorem 13.91 ob- 
taining an optimal control (c* ,ir* , a*) for the original problem. The result again follows 
from arguments exploiting the Markovian property of our controlled system, we only state 
it without giving a proof for sake of brevity. 

Theorem 3.10 Assume that the hypotheses of Theorem 1 3. iA hold. Then an optimal control 
(a* , c* , 7T*) is given in closed loop form as 

a* k G argmax ^ R V(0,Rr k - a,a), k > 0, 



<c* s = C*(s-r k ,X. 



< = n*(,-rar^- a -,^), 



S G [T k ,T k+1 ), k>0, 



s G [T k ,T k+ i), k>0. 
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When U is strictly concave we have uniqueness of optimal controls. 

Proposition 3.11 Assume that the assumptions of Theorem \3.9\ hold. If U is strictly 
concave, then V and V(t, •, •) are strictly increasing (V with respect to both x,y) and strictly 
concave, and we have uniqueness of optimal controls. 

Proof. The fact that U is strictly concave implies also that it is strictly increasing. This 
yields also strict monotonicity of the value function V (using a standard argument that can 
be found, e.g., in [26j Prop. 2.1] or [61 Prop. 4.7]). 

Now, let < 7*i < r2 and take two optimal controls (c , it , a ), (c 2 ,7r 2 ,a 2 ) for r%,r2 
respectively, their existence being provided by Theorem l3.101 Since V is strictly increasing, 
it must be c 1 ^ c 2 . Let A G (0, 1) and set r\ = Xr\ + (1 — A)r2 and 

(c x ,ir x ,a x ) = A(c 1 ,7r 1 ,a 1 ) + (1 — A)(c 2 , 7r 2 , a 2 ). 

Then by linearity of the state equation ©-([U) we have (c A ,7r A ,a A ) E A(r\). On the other 
hand, since U is strictly concave and c 1 ^ c 2 , we get V(r\) > XV(r\) + (1 — X)V(r2), which 
proves that V is strictly concave. 

Strict concavity of V implies also strict concavity of G 7 [y](i,-,-) for every t > 0. In 
the same way strict monotonicity of V implies strict monotonicity of G 7 [V] with respect to 
x, y, which implies strict monotonicity of V with respect to x, y. Then, arguing as before, 
one gets strict concavity of V(t, •, •) too. 

Let now t,x,y > and suppose that (c 1 ,^ 1 ), (c 2 , 7r 2 ) are optimal controls for the 
auxiliary problem starting from (t,x,y). As a consequence of strict concavity of U, we get 
c\ = c\ =■ c*. As a consequence of strict concavity of G 7 [y], we get X t ' x ' c *' 7r * = X*> x,c * ,7r 2. 
From that, one easily derives ir\ = tt^. 

Now consider the original problem with initial value r and take two optimal controls 
(c ,7r , a ), (c 2 ,7r 2 ,a 2 ). Strict concavity of V implies that a i— > V(0,r — a, a) has a unique 
maximizer in [0, r], so ckq = a§ =: Qo- On the other hand the controls (c , 7r )|[o )ri ) and 
(c 2 , vr 2 )|[ n \ are (locally) optimal for the auxiliary problem in [0,ti) with initial value 
(x,y) = (r — ao,ao). By uniqueness of optimal controls for the auxiliary problem, it must 
be (c 1 ,tt 1 ) = (c 2 ,7r 2 ) in [0, ri). By the Markov property, the same argument applies to 
generic intervals [Tfc,Tfc+i), so the proof is complete. □ 

4 Power utility 

In this section we consider the problem when the utility function is 

U(c) = -, P e (0,1). 
p 

In this case the Legendre transform of U is the function 

1 — p p 

U(w) = w 1 -p , w > 0. 

V 

We notice that in the case p < the problem is investigated in [I], but supposing full 
observation of the illiquid asset, corresponding to the case 7 = 1 in our setting. 



19 



4.1 Reduction to one space variable 

Due to the fact that the utility function U is homogeneous, the function V is a power 
function and the function V is homogeneous in the variables x,y. So (|38p can be reduced 
to a PDE involving just a one-dimensional state variable. Indeed, since 

A(£r) = £A(r), V£ > 0, Vr > 0, 

taking into account the fact that U is homogeneous of degree p, it is straightforward to 
show that for some $o > 

V{r) = <Z> r p , Vr > 0, (52) 

Moreover, since we have also 

At(fr) = ZMx), V£ > 0, \fx > 0, 
taking into account (|52p and (|93|) . it is straightforward to show that 

V(t, £x, £y) = e V{t, x, y), V£ > 0, Vx, y > 0. (53) 
Hence, following e.g. [23], we perform the change of variables^ 

z = -, x>0, y>0. (54) 

By this change of variable and (|53|) , we can rewrite V in separated form as 

V(t,x,y) = y p <P(t,z), t>0, z>0, y > 0, (55) 

where 

$(t,z) = V(t,z,l), t>0, z>0. (56) 



Set, for y > 0, 



~ (57) 

y y 



Denote by A^ p (IRj_,M) the space of measurable functions ip such that \ip(t, z)\ < C(l + \z\) p 
and consider the nonlinear functional 

n : M P (M. 2 + ,R) — )■ R, 

4,(0, z) 



ip i — > Ho[ip] ■= sup 



z > (1 + z)p 



Then one can check that V(r) = H[V](r) = H [$]r p , so $° = H [$] and setting 

p(t, Z ) = cr[t^e\(t,z,i) 



s Here we restrict the analysis to the case, y > 0. In the case y = 0, the problem is already a one- 
dimensional Merton type problem for which one can deal with explicit solutions to the HJB equation (see 
Subsection 14. 4[) . A necessary and sufficient condition ensuring that y — is never optimal is given in 
Proposition 14.81 below. 
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and taking also into account (|93p we get 

G^[HV](t,x,y) = y p G"<lHV](t,z,l) = f (t, z)Wo[*]. (58) 

Plugging (|53|) -- (|58|) (with the formal argument ip in place of $) into (|38|) and dividing by 
y p , we get an equation for $ in M+ x (0, +00): 



tp t + (/3 + A)y> - \p(t,z)H [ip\ - sup H^z,?,^,^) 

c>o, seR 



0. 



(59) 



where 



Pz,Vzz; c, 7r) = (p 2 +7 2 (!-p 2 ))y {p{p-l)^>+2(l-p)zip z +z 2 ip z ^ +b Y (pp-z(p z ) 
+ U(c)- ap z + 7r ((&£ + pdia L {p - l))<p x - pcria L ZLp zz ) + Tr 2 -^(p zz , 

We notice that "Hofv] is a nonlocal term in (159p . 

By monotonicity and concavity of V(i, •> •) we deduce also monotonicity and concavity 
of $(t, •). Hence, we have & z > 0, $ Z2 < in the sense of sub(super)differentials. Now we 
introduce a change of variables to simplify the computations, denoting 

a ~ pai 

= TT Z, 

so (I59j) can be rewritten after suitable simplifications as 

K 2 

- <p t + K\tp - K 3 z(p z - Xp(t,z)Ho[<f] 7T Z Vzz ~ SU P Hl v (ip z ,(p zz ;c,9) = 0, (60) 

1 c>0, 



where 



and 



1 



U{c) - ap z + 9K^ Z + ^&K%<p 



(61) 



K A = /? + A + — ^-^(1 - p) - 



0"L 



7 2 P fe / 



Ki = b L - paia L (\ -p), 
K 2 = cr L , 

K 3 = 7 2 f-6 / + ^ + (l-p 2 )(l-pW 



(62) 



^4 



-any/1 - p 2 



Since & z > and <3? zz < in the sense of sub(super)differentials, we can rewrite ()60|) as 

-^ + ^-if 3 ^-Ar(t,^oM-^) + ^— -?$z 2 <Pzz = 0, (63) 

2 A| ip zz 2 

with the convention 



h 




0, if = 0, 
+00, if |/i| > 0. 
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By (f35j) we get the boundary condition 

POO f-OO 

$(t,0) = Ho[$] / e- K ^ s -Q\p(s,0)ds = $° / e- K ^ s - t) Xp(s,0)ds, (64) 

Jt 

and by (|37p we get the growth condition 

$(t,z) < Kye kj <v\l + z) p . (65) 

Due to the results of Section 3 and to the argument above, we get the following (according 
to the definition of viscosity solution of [1] ) : 

Proposition 4.1 The function <3? is continuous onM.^_, strictly increasing and strictly con- 
cave in z, and is the unique viscosity solution over R+ x (0, +oo) of (|6Up fulfilling the 
boundary condition (f64"|) and the growth condition ([65]) . 

4.2 Smoothness of the value function 

In this subsection we show that the value function is smooth. We can freeze the nonlocal 
term in (|59p and, by using Proposition 14.11 and standard comparison results for viscosity 
solutions of second order parabolic PDE's (see, e.g., [3]), we get the following: 

Proposition 4.2 <I> is the unique viscosity solution over M + x (0, +oo) of 

K 2 

- ft + K\f ~ K 3 zip z - A/ 7 (t,z)$° - -^-z 2 Lp zz - sup H cv (z,ip,<p z ,(p zz ) = 0, (66) 

z £>o, 0eR 

fulfilling the boundary condition (|64l) and the growth condition ([65]). 

The nature of the HJB equation (|66p is sensitive to the value of 7. More precisely one has 
to distinguish, from the point of view of PDE's theory, the cases 7 = and 7 7^ 0. Indeed, 
in the case 7 7^ we have K4 > and the PDE is a fully nonlinear nondegenerate parabolic 
equation, while in the case 7 = we have = and the PDE is degenerate. While a good 
regularity theory of solutions is available in the nondegenerate case, in the degenerate case 
only the theory of viscosity solution applies in general. In the latter case, the possibility 
of getting good regularity results (like classical solutions) strongly relies in the specific 
structure of the equation. In both cases we can prove that the solution is smooth enough 
to apply Theorem 13.91 Due to the considerations above, the cases 7 = and 7 7^ must be 
treated in a different way: the degenerate case 7 = requires a passage to a dual problem, 
which is explained in detail in [8]; on the other hand in the nondegenerate case 7 / 
the known regularity theory can be used, but needs some technical nonstandard results to 
localize the equation and restrict the set of controls to a compact one. These results, on 
which relies the proof of the following regularity result, are provided in the Appendix. 

Theorem 4.3 $ G C7 1 ' 3 (R+ x (0,+oo);R) and <5> z G C 1 ' 2 (R+ x (0,+oo);R). Moreover 
> 0, $ zz < in R + x (0, +00). 
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Proof. As we said, we prove the claim in the case 7^0, referring to [8] for the case 7 = 0. 
We note that in the latter case we get more reg ularity: indeed $ G C 1,00 (M+ x (0, +00); R) 
and all the space derivatives lie in the same space. 

Given (i,z) G R+ x (0, +00) and e G (0,z), consider the set D e (i,z) defined in (jlOip 
and let V(D £ (t, z)) be the parabolic boundary of D £ (t, z) defined as 

V(D £ (i,z)) := {i + e} x [z - e,z + e] U [i - e,i + e] x {z - e,z + e}. 

Then, by Proposition IA.6I and |10|, Ch.V, Cor. 8]), is the unique continuous viscosity 
solution on V(D £ (i, z)) of the HJB equation (|112|) - which is the same as (|66|) but with 
constraints on the set of the variables c, 9 - with Dirichlet continuous boundary condition 

<p = on V(D £ (t,z)), (67) 

On the other hand, by p21 Th. 3, Sec. 6.4, p. 301] (see also Example 8, Section 6.1, p. 279, 
of the same book) there exists a solution C 1,2 (D e (i, z); K) of (|112|) with boundary condition 
(|67p . Since such solution must be also a viscosity solution, we conclude that 



$ G C 1,2 (.D e (M);R). (68) 

Moreover, by Lemma IA.5I we also have for (t, z) S D e (t,z) 

*x(t,z) > m e > 0, *«(t,«) < -4 < 0. (69) 

Due to (|68|) - (|69|) and Proposition 14. li and by arbitrariness of (i, z) S R+ x (0,+oo), we 
get that $ G C 1 ' 2 (R+ x (0, +oo);R) is a classical solution of (j63[) and $ z > 0, < in 
R+ x (0,+oo). 

Now we apply Lemma I A. 71 which allows us to formally differentiate (|63p and conclude 
that $ z is a viscosity solution tc0 



5t + (K X + K 3 + ^g+(K 3 - Kl)zg z + U'(g)g 
K 2 K 2 n 2 \ 

^z 2 + ^ 9 -)g zz + \^f z = 0, (70) 

with Dirichlet continuous boundary condition 

g = $ z , on V(D £ (i,z)). (71) 

Again by |101 Ch.V, Cor. 8.1], the function & z is the unique viscosity solution to this 
problem. On the other hand equation (I70p is a quasilinear uniformly parabolic equa- 
tion, so it admits a solution of class C l,2 {D £ (t, z); R) (see, e.g., [18} Th. 12.22] - with the 
assumptions of Th. 12.16 of the same book). As above, by uniqueness we deduce that 
<fr z G C 1,2 (D £ (t, z);R), hence in particular that $ G C l ^(D £ (t, z)\ R). By arbitrariness of 
(t, z) G R+ x (0, +00), we get the final claim. □ 



9 Although U'(g)g z and g 2 /gl are not well-defined for g = or g z = 0, we may use (|69p to replace these 
terms by bounded continuous functions of {g,g z ) coinciding with them whenever m e < g < M e , g z < —S E . 
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4.3 An auxiliary closed loop equation 

In this subsection we study a closed loop equation associated to the feedback maps provided 
by the maximization in the HJB equation (fSUj) (for z > 0), and by the fact that ^(0) = 
{(0,0)} (for z = 0): 



[0, ifz = 0, [0, if 2 = 0. 

Of course these maps are measurable in the couple (s, z). Moreover, due to Theorem l4.31 the 
maps C*(s, •), 0*(s, ■) are locally Lipschitz continuous in (0, +oo), uniformly in s £ [0,T] 
for all T > 0. We associate to these maps the closed loop equation 

dZ s = -C*(s, Z s )ds + 6*(s, Z s ) (fcxds + K 2 dW s S j + Z s (K 3 ds + K^dB^) ^ 
2* = ^, 

where K x = b L - poio L and ^ 3 = 7 2 (-6/ + + (1 - p 2 )af}. We note that Ki and 

are not the constants appearing in the HJB equation (|60p . We comment on that in 
Remark 14,51 below. 

Proposition 4.4 Given (t,z) € there exists a unique solution Z*' z '* > to (1721. 

Proof. Existence. If z = the claim is clear, just by taking Z t,z '* = 0. Let z > and 
T > 0. Due to local Lipschitz continuity properties of C*(s, ■), 0*(s, •), using standard 
SDE's theory (see, e.g., [TBI Ch. 5, Th. 2.9]), we get for each e G (0, z) the existence of a 
unique solution Z t)Z)S € [e, e" 1 ] in the stochastic interval [i, rj), where is implictely 
defined in terms of the solution itself as 



rj = inf {s e [t,T] \ ZY' S < e or Z f s ' z ' £ > e" 1 }, 
with the convention inf = T. Of course, if e < e', we have > and 

z t,z,e _ z m, £ ' on ^ r T^ V0<e<e'. (73) 

Set 

r T = liniT T . 

e4.o 

Then by f|T3[) there exists a unique solution Z*' 2 '* > in the interval [t,r T ). We now show 
that this solution can be extended to the whole interval [t, T], which, due to the arbitrariness 
of T, will imply the claim. By a Girsanov transformation, there exists a probability Q T 
locally equivalent to P, and Q T -Brownian motions W, such that (|72p may be rewritten 

as 

dZ s = -C*(s,Z s )ds + K 2 e*(s,Z s )dW s + K 4 Z s dBp. 

By nonnegativity of C* and Z t,z '*, the process Z t,z '* is a nonnegative Q T -supermartingale 
on [t, t t ). It can be extended to a Q T -supermartingale (L 1 bounded) on [t,T] by setting 
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Z t,z '* = in [t t ,T]. Hence, by Doob's convergence Theorem (the usual proof for deter- 
ministic intervals - see e.g. Theorem 6.18 in [15] - can be easily adapted to our stochastic 
interval [t,r T )), there exists a finite random variable Z*'*'* such that 



lim Z\> z >* = Z z f , Q J -a.s.. (74) 

S /-T T T 



Since Q T ~ P, we also have 



lim Z*' 2 '* = Z 1 ' 2 '*, P-a.s.. (75) 

S /-T T 

Immediately (|75p yields the desired extension on {t t = T}. Let us now consider the set 
{t t < T}. On this set we have Z*'^'* G {e, e -1 }, so that by ([75]) necessarily Z*'^'* = 
almost surely therein, getting 

lim Z*' 2 '* = a.s. on {r r < T}. (76) 

Therefore, we may now extend Z t,z '* to a solution defined over [t, T] on {r T < T} by setting 

Z*' 2 '* = 0, for s G [r T ,T]. 

Uniqueness. The solution is clearly unique on the stochastic interval [t, t t ) defined in 
the existence part. On the set {t t < T}, when it reaches 0, it must stay there, since 
it is a nonnegative Q T -supermartingale. Therefore we have uniqueness on [t, T] and, by 
arbitrariness of T, on [t, +oo). □ 

Remark 4.5 The constants K\ and that we take in (|72|) are different from the constants 
K\ and K% in the HJB equation (|6U|) . These are indeed the right constants that allow to 
come back to the auxiliary control problem in the couple (X,Y) of Section 3, as we will 
show in the next subsection. 

To explain better this fact, we note the following: 

(1) What we did in the previous subsection was to reduce the solution to the HJB in two 
variables to the solution of the HJB in one variable, using homogeneity of U and 
simplifying the variable y. 

(2) It is possible to do this reduction directly on the control problem. To do this one has to 
define the new state process Z as the ratio of the old state variables, find an equation 
for it and rewrite the utility functional in terms of it. Concerning the equation for 
Z, let x > 0, y > 0, let (c,vr) G A (x) and set X = X^ x ^ , Y = Y ^. Denoting 
Z s = one may check that Z solves the SDE 

dZ s = -c s ds + 6 s (K 1 ds + K 2 dW s ) + Z s (K 3 ds + K 4 dBP), (77) 

where c s = and 6 S = ^ — ^-Z s , so with the same constants as the ones appearing 
in (|72p . Then, by homogeneity of U one may rewrite the objective functional of the 
auxiliary control problem of Section 3 as 



E 



' />oo 

/ e-^+» s (Y s y (U(c s ) + \(P\y](8, Z s , 1)) ds 
Jo 



(78) 
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Similarly to what we did for the HJB equation, we now need to "simplify" the term Y p . 
The natural way to do that in this context is to change probability. Indeed consider 
the probability P corresponding to the density process E ^ ) P | • Under this probability 
the processes 

W s := W s - ppa lS , BP := B« - - p 2 s 
are standard Brownian motions, (|77p writes equivalently as 

dZ s = -c s ds + e s (K 1 ds + K 2 dW s ) + Z s (K 3 ds + K i dBi 1) ). (79) 
and (1781) can be rewritten as 



y p E* 



oo 

e- K " s (U(c s ) + Xf(s,Z s )ds) 



o 



(80) 



As expected, the HJB equation associated to the control problem with state equation ([79 
and objective functional (j80j) is exactly (|60|) . 



4.4 Back to the original problem 

As a consequence of Theorem 14.31 and Proposition 14.41 we can apply Theorem 13.91 Indeed 
we have the following: 

Proposition 4.6 IfU(c) = c p /p, all the assumptions of Theorem \3.9\ are satisfied. 

Proof. (|43p (i) follows from the smoothness of <3? (Theorem 14. 3j) and (|55p . Concerning 
(1431) (ii) . we observe that, when y = 0, we get from (f30l) . (|52j) and (i24l) 



V(t,x,0) = sup E 



(81) 



Since this is a standard homogeneous Merton type problem we get V(t, x, 0) = Kqx p for 
some Kq > 0, which proves in particular (03]) -(ii). 

Let us show now that also the assumptions regarding the existence of the feedback maps 
and the solvability of the closed loop equation in Theorem 14.31 are satisfied in this case. Set 



x\ pa! 

n = 9 ( s, - H z, 

y' o- L 



and, for any x > 0, y > 0, 



C*(s,x,y) = yC*(s,-), U*(s,x,y) = yU*(s,- 
\ yJ \ y 

The maps C*,H* satisfy the assumptions (|44p -(|45 p of Theorem 13.91 Moreover, one can 
check by integration by parts that Z is a solution to the closed-loop equation (j72|) if and 
only if X := ZY t,y is a solution to the SDE (|46p . So, also the assumption of Theorem 
3.91 about existence and uniqueness of nonnegative solution to the closed loop equation is 
satisfied due to Proposition 14.41 The case y = can be treated separately, and more easily, 
working on the degenerate Merton type control problem (|8ip . □ 
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4.5 Two properties of the optimal allocation in the illiquid asset 

Here we prove two important properties of the optimal allocation strategy in the illiquid 

x* 

asset. The first one concerns the fact that the ratio y^r- is constant. 

T k 

x* 

Proposition 4.7 If U(c) = c p /p, the optimal rebalancing proportion y^r- at the trading 

T k 

times of the illiquid asset does not depend on the current value of the wealth. Indeed 

K _* . m*) 



where the value z* above is well defined under the convention that z* = +00 if the supremum 
°f (i+z)p ^ s n °t attained, in which case there is no investment in the illiquid asset. 

Proof. First of all, we observe that, by strict convexity of $ the function a 1— > a p $(0, ^ — 1) 
has a unique maximum point on [0, r], which we call a*(r). Then, due to the homogeneity 
of V, we have 

X *k R r k ~ a *( R r k ) 

Y* k a*(R Tk ) ■ { ] 

Now we notice that by definition of a* 

R Tk 1 . *(0,^) 



1 = argmax 2>0 



a*(R Tk ) ° - {l + zY' 

which yields the claim. □ 

The second property is concerning with the possibility of having nonzero optimal investment 
in the liquid and illiquid asset. 

Proposition 4.8 

1. olq < r (if and only ifr>0). 

2. a* >0 if and only if £ > ^0 

Proof. 1. Of course if r = then one needs to have ag = due to the state constraint. 
Instead assume that r > and assume, by contradiction, that Oq = r. Then this would 
yield z* = 0, hence, due to Proposition 14.71 above, a* k = R Tk for all k € N. We should 
conclude that the liquid asset would be always 0, implying by the state constraint c* t = 0, 
which cannot be optimal as V(r) > 0. 

2. Necessity. Consider the Merton problem described in Remark 12.61 where the agent 
invests continuously in L and /, with the constraint that the proportion invested in / be 
in [0, 1], and with value function denoted by V M ' 2 . The optimal investment proportions in 
L and / for this problem are given by : 

(u* L ,u}) = argmax \p(u L b L + uj6j) - — ^ {u\a\ + ujaj + 2pu L uia L ai) 1 
«l eR, ?ire [0,1] I 1 J 



10 This condition is the same as the one in the Merton (liquid) problem with two assets (the same result 
is obtained in [I] in the case of full observation). 
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. In 



Taking first the supremum on Hi, one can see that uj = if and only if bj < 
this case, denoting by V A1,1 the value function for an agent investing only in L, we have 
V M,2 = yM,i g ince obviously V M ' X < V < V M ' 2 , we obtain V = V M '\ and the optimal 
strategy for our original problem never invests in the illiquid asset I. 



Sufficiency. Assume — > By the homogeneity property due to the power utility, 



it is enough to show that h'(0 + ) > 0, where h(a) = 7(0, 1 — a, a). Noting that 



h'(0 + ^ 



lim 

7}->0 

lim 



7(0,1-^^-7(0,1,0) 



(1-7/)* 



r?->0 rj 



7(0,1, 



1 — 7]' 



7(0,1 + 



1 — T] 



0) 



lim (1 - r/f- 1 



lim -(7(0,1,5) -7(0,1 + 5,0)) 



7,(0, 1,0+)- 7,(0, 1+0), 



we will show that the latter is strictly positive. 

Consider the auxiliary problem with initial data (t,x,y) = (0,1,0). In this case the 
problem is the Merton type problem (|8ip with homogeneous value function 7(0, x, 0) = 
Kqx p , so 7 x (0, l + ,0) = pKq. By solving the HJB equation for this problem, one can see 
that Kq is the unique positive solution to 



/3 + A 



pb\ 



2(1 -PK 



Kq - (l-p)p *-pK q 



V 

l-p 



A$ c 



(83) 



and the corresponding optimal wealth process X* is given by 



dX* 



-Cfdt + 7T 



Lt ' 



(84) 



where 



4 = P -—K Q x: 



t ■ 



7T/ 



b L 



*i(i-p) 



X*. 



(85) 



Considering an agent with initial wealth (1,5) following the same investment /consumption 
strategy we get 



7(0,1,5) > E 



-(/3+\)t 



(U(c* t ) + \G[V}(t,X?,Y t °' S ))dt 



L J0 



Therefore 

9(0,1,5) -7(0,1,0) 
5 



> *E 


" POO 








Jo 




poo 


= A<D° 


L 







(G[V](t,X*,Y t ' 6 ) - G[V](t,x;,0))dt 



-(/3+A)t E 



(xt + Y^j t y-(x;r 

5 



dt. 
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Letting 5 — > 0, applying Fatou's lemma, and observing that Y°' S Jt = 5It, from the inequal- 
ity above we get 



v y (o,x,o + ) > P <s>°\ / e - (/3+A)4 E {x;y- l i t dt 

Jo L 



= P ^ x r ew r {x ^_ {bl _p^i ))t ) dt 

J0 V K <*L ' 

> pK = K(0,l+,0), 

where in the equality in the middle we have used ([T]), ([2]), ([83]) . ([M)) and ([85]) . and in the 
strict inequality we have used ^f- > The proof is complete. □ 

The result above says in particular that, when the two assets are uncorrelated, there is 
investment in the illiquid asset even if the Sharpe ratio of the liquid asset is higher of the 
one of the illiquid asset, i.e. there is still diversification in the allocation of the portfolio 
between the two assets. 



4.6 Numerical approximations 

In this subsection we present an iterative scheme to approximate <I> and <£. This procedure 
is illustrated more extensively in the case 7 = in [11] , where the value functions V and 
V are approximated. Here we describe such procedure for 7 £ [0, 1] with regard to ( I >0 and 

i.e. for the power utility case, as it will be used in the next section to produce numerical 
results and comments in this special - but relevant - case. 

Let us look at (159|) . Because of the nonlocal term T~Lq [ip] , we cannot approximate directly 
the value function $ as viscosity solution of a PDE, but we need to define an iterative 
scheme. Fix T > 0. Starting with 

<D°>°< T = 0. (86) 
we define inductively the sequence ($°' n,T , $ n ' T ) as follows: 

- Given n > and 3>°> n '' r 5 we define <& n ' T on as the unique (constrained viscosity) 
solution on [0, T] x IR+ to 

- + {p + A)$"' T - \^ T f {t, z) - sup H cv (z, <D"> T , $ n z > T , $ n J; c, tt) = 0, (87) 

c>o,#eR 

with boundary condition 

$ n > T (t,0) = ^ T j\- K ^ s - t) Xf'(s,0)ds,, (88) 
and terminal condition 

<S> n ' T (T,z) = 0. (89) 

- Given n > and <]? n > T , we define $ > n+1 > T by 
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Then one can prove (see [TT] for details) that 

lim ^+l,T = $0) Hm $ n+l,T = $ 

n—too, T— >oo n— >oo, T— >oo 

The rate of convergence above is sensitive to the value of A (see [7]): 

- The bigger is A, the slower is the convergence in n. 

- The smaller is A, the lower is the convergence in T. 

5 Discussion 

In this section we provide and discuss some numerical experiments we performed in the case 
of power utility by means of the iterative approximation procedure that we have described 
in Subsection Many more tests and more discussion could be done looking deeply at the 
properties of the optimal paths0 Since the purpose of this paper is mainly methodological, 
here we limit ourselves to discuss some key features to show how our methodology can be 
applied. 

We have taken as values of the parameters 

/3 = 0.2, p = 0.5, b L = 0.15, a L = 1, 6/ = 0.2, oi = 1. 

Let us explain our choice for these parameters. First, the illiquid asset should have a higher 
Sharpe ratio than the liquid one. This is economically intuitive, and moreover ensures that 
for any value of the correlation p, it will always be optimal to invest something in the 
illiquid asset. Second, we want that the optimal investment proportion in / be in [0, 1], so 
as to observe mainly the impact of the constraint on trading dates rather than the impact 
of the proportion constraint (induced by discretization). As in \21\ I12j this is obtained for 
parameters corresponding to a high risk-return market. 

We then vary the other relevant parameters 7, A, p, representing respectively the ob- 
servation of the illiquid asset between two trading times, the liquidity of the market, and 
the correlation between the two assets. We have solved the PDE's (|87p using an explicit 
finite-difference scheme for parabolic viscosity solutions (see [10^ Ch. IX]). More precisely, 
this has been done after a change of variable 

M+ -> [0,1), z k> z = — = , 

2+1 x + y 

inducing a corresponding transformation $ 1— )■ <&, in order to work on the bounded domain 
[0, 1). We have taken T between 1 and 5 (depending on A according to what said at the end 
of Subsection l4.6p and used a uniform grid on [0, T] x [0, 1] with time step length 5- 10 -4 and 
space step length 0.02. The numbers / 7 (t,z) were computed beforehand at each point of 
the grid using an L 2 -optimal quantization grid for the gaussian law with N = 5000 points. 
Finally, the derivatives have been approximated by finite difference. 

See e.g. the tests performed in [T]. 
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5.1 Value function and cost of illiquidity 

In this subsection we study the cost of illiquidity by looking at some quantities related to 
the value function V. 

The first (absolute) way to measure the cost of illiquidity is, of course, by looking at 
the differences of the value functions corresponding to different values of A. In Figure Q] 
we represent the value function Wlfin as function of the correlation parameter p and with 
fixed 7 = (no observation between trading times). The different lines correspond to 
different values of the liquidity parameter A. We also sketched the graph of the constrained 
and unconstrained Merton problem with the two assets L, /, with / considered as liquid 
as well. We observe that the illiquidity, measured by A, has a considerable impact on 
the value function for negative correlation. This is expected, since to deal with negative 
correlation the agent is naturally driven to hedge the fluctuations of one asset by assuming 
an opposite position in the other one. However, this is not always possible when one has 
constraint on the strategies. This effect can be already seen comparing the constrained and 
the unconstrained Merton model, and it is amplified by the presence of illiquidity, which 
induces further constraints on the strategies. To this regard it is worth to stress that we 
cannot expect to have always convergence to Merton unconstrained problem for A — > oo. 
Indeed, the presence of illiquidity, even in the case of high A - high trading frequency - 
immediately induces a constraint on the investment strategies in I in order to satisfy the 
state constraint - see (J7|) . This fact may produce a gap between the fully liquid case and the 
illiquid case (even if the trading frequency is very high) . Actually, the limit case for A — > oo 
corresponds to the constrained fully liquid Merton problem, i.e. the Merton problem with 
the constraint that the investment a = (ctt)t>0 m I does not admit borrowing or short 
selling: < at < Rt for every 

On the other hand, the impact of the observation parameter 7 on the value function 
V(l) is observed to be low, both in absolute value and in percentage change. This low 
impact can be observed without regard to the liquidity parameter A and to the correlation 
parameter p, as shown by Table Q] for the extreme cases 7 = 0, 1 for p = 0. We performed 
the same analysis for other risk aversion (not reported here for brevity; we did it also in 
the case of higher risk aversion, in particular for negative values of p) as well as for other 
values of the correlation p, and still we have observed a low impact of 7 on the results. 
We observe that the change of 7 has naturally negligible impact both when A — > and 
A — > +00. Indeed, when A = the illiquid asset is not traded; when A = +00 the time 
distance between trading times goes to 0, so the observation between trading times plays 
no role. So, as we expect, the impact of 7 no the value function is higher for intermediate 
values of A. 

For comparison purposes we have also computed the value of V(l) in the model of |21j . 
i.e. with no liquid asset and no observation on the illiquid asset (and with the same values 
of the parameters hi and 07). We see that in this case adding the illiquid asset, V(l) grows 

12 For another initial wealth r the same measure can be performed by the relation V(r) = V(l)r p . 

13 Of course, when the optimal solution of unconstrained Merton problem satisfies this constraint, con- 
strained and unconstrained Merton problems are equivalent and we have the convergence to the (uncon- 
strained) Merton problem for A — > 00. 
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1.6 

-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



A = 1 A = 50 

A = 5 constr. Merton 

A = 10 unconstr. Merton 

Figure 1: Value function V(l) as a function of p with fixed 7 = 0. 



A 


1 5 10 50 


Merton 


7 = 


1.66641 1.70493 1.71257 1.71945 


1.72133 


7 = 1 


1.66995 1.71121 1.71656 1.72036 


1.72133 



Table 1: V(l) for various 7, A and fixed p = 0. 



of about 3% without relevant changes with respect to A. 



A 


1 5 10 50 


Merton 


V(l) 


1.61973 1.65377 1.65987 1.66526 


1.66667 



Table 2: V(l) for various A in the case of no liquid asset (and 7 = 0). 



Another way of measuring the cost of illiquidity is to define it - see [21] - as the extra 
amount of initial wealth e(r) needed to reach the same level of expected utility as an investor 
without trading restrictions and initial capital x. Hence, it is then computed as the solution 
to the equation 

V(r + e(r)) = V M (r), 

where Vm is the value function of the corresponding unconstrained Merton problem. 

Tables HUE] reproduce the value of e(l]H| for different values of A and for 7 = 0, 1, 
respectively for p = 0, p = 0.5, p = —0.5. As in the case of 1^(1), the impact of the 
parameter 7 is low in an absolute value, reaching its maximum for intermediate values of 
A. Moreover, we also observe that the relative impact — — e*=i(i) — ~' *' e ' ^ e P ercen tage 

14 In our setting of power utility, the cost of liquidity e(x) is proportional to x. We therefore study the 
cost of liquidity per unit of initial wealth r = 1. 
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change of e(l) passing from 7 = 1 to 7 = 0, is much higher than in the case of V(l). 



A 


1 5 10 50 


7 = 


0.067 0.0193 0.0103 0.00218 


7 = 1 


0.062 0.0119 0.0056 0.00112 



Table 3: e(l) for various 7, A and fixed p = 0. 



A 


1 5 10 50 


7 = 


0.0337 0.00892 0.00462 0.00095 


7 = 1 


0.0303 0.00491 0.00237 0.00051 



Table 4: e(l) for various 7, A and fixed p = 0.5. 



A 


1 5 10 50 


7 = 


0.2511 0.1127 0.0700 0.0161 


7 = 1 


0.2493 0.1030 0.0614 0.0120 



Table 5: e(l) for various 7, A and fixed p — —0.5. 



Also in this case we have performed the analysis on the impact of the parameter 7 on 
the results for other for other risk aversion, observing a low impact. 

5.2 Optimal policies 

Now we look at the optimal policies. Hereafter, without loss of generality (as we are in the 
power utility case) we assume that the initial wealth is r = 1 - hence, investment expressed 
as amount coincides with investment expressed as portfolio proportion at time t = 0. 

In Figure [2] below we represent the optimal allocation in the illiquid asset as proportion 
of the wealth - i.e. the quantity z = 1 — z* - as function of the correlation p for 7 = (no 
observation between trading times). The different lines correspond to different values of A. 
When A is very low, clearly the optimal investment proportion in the illiquid asset is close to 
0; on the other hand, when A is very high, it is close to the one of the constrained Merton 
problem. In this last case the investment in the "illiquid" asset is higher for higher (in 
absolute value) correlation, as expected. So, for increasing values of A, the corresponding 
graphs lie between these two extreme cases, clearly increasing with A. in particular, the 
optimal investment in the illiquid asset will be lower for lower liquidity, which is quite 
intuitive. 

The impact of the observation parameter on the optimal allocation proportion in the 
illiquid asset z* is represented in Table [6] for different values of A considering the extreme 
cases 7 = and 7=1 and setting p = 0. Clearly the agent will invest more in the illiquid 
asset if he can observe it continuously. Moreover the impact of 7 is negligible in the extreme 
cases A = 1, A = 50 and of the order of 6% when A = 3, 5, 10. So, as for V{\) and e(l), the 
impact of the parameter 7 is higher for intermediate values of A. We performed the same 
analysis for other values of p and observed similar results. 
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-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 

A = 1 A = 50 

A = 5 constr. Merton 

A = 10 unconstr. Merton 

Figure 2: Optimal investment proportion z in the illiquid asset as function of p for 7 = 0. 





A = 1 A = 3 A = 5 A = 10 A = 50 


Constr. /Unconstr. Merton 


7 = 
7=1 


0.18 0.3 0.34 0.36 0.4 
0.18 0.32 0.36 0.38 0.4 


0.4 
0.4 



Table 6: Optimal investment (proportion) in the illiquid asset z* with p = 0. 



Let us go now to the analysis of the feedback maps defining the optimal consumption 
and the optimal investment in the liquid asset. We want to see the influence of the portfolio 
proportion invested in the illiquid asset on these maps, so it is meaningful to study them 
with respect to the variable z = Hence, we look at the functions 

C*(t,z) := zC* (t, J — 1 J , U*(t,z) := zft (t, ~ - l^j . 

indeed, by homogeneity, given the total amount x + y, we have 

C*(t,x,y) := (x + y) C*(t,z), U*(t,z) := (x + y) fi*(t, z). 

So if we set the total initial capital equal to 1, the functions C*(0, z), II* (0, z) yield re- 
spectively the consumption and the optimal amount invested in the liquid asset at time 
0. We therefore plot the graphs of C*(0, z), II*(0, Note that we must pay attention 
in comparing these maps for different values of the parameters 7, p for t > 0. Indeed they 

1 In our model, in practice at time t = the agent following an optimal strategy will always have fixed 
proportion £* invested in the illiquid asset, as we are assuming t — being a trading date. Anyway, it is 
still meaningful to look at the graphs of C*(0, •), II* (0, •), as this corresponds to suppose that the t = is 
not a trading date - see [T][24] and Remark l3.3f iii) - and one can be interested in the analysis of the optimal 
consumption and of the optimal investment in the illiquid asset as function of the initial repartition z. 
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are defined on the variable z, which refers to the stochastic process Z = 1/(Z + 1). But Z 
itself depends on the parameters 7, p, so the feedbacks maps C*,IL* do not read the same 
input for different values of the parameters 7, p. Nevertheless, at time t = the meaning 
of this variable is clear, as it is then simply the proportion of the total wealth invested in 
the illiquid asset. Computing these functions for different values of 7 we observe very small 
changes in the shape and in the values (less than 5 % relative change for all values), so we 
plot their graphs only in the case 7 = 0. 

0.4 |— 
0.35 

0.3 - 

0.25 - 

0.2 - 

0.15 — 

0.1 — 

0.05 — 

— 

0.2 0.4 0.6 0.8 1 

A = 1 A = 10 Merton 

A = 5 A = 50 




Figure 3: Optimal consumption rate C*(0, •) as function of z with p = 0, 7 = 0. 

In Figure [3] we plot C*(0, •) only for p = 0, as different values give similar shapes. As 
in [12] we observe that the influence of A on the optimal consumption rate depends on 
the proportion of illiquid investment z: when z is close to 1, hence most of the portfolio 
is constituted of illiquid wealth, the investor faces the risk of "having nothing more to 
consume". As a consequence, the further away the next trading date is, the smaller the 
consumption rate should be, thus C* is increasing in A. When z is far from 1, the opposite 
happens: when A is smaller, the investor will not be able to invest optimally to maximize 
future income and should consume more quickly. 

In Figures U]-[5]-[6j we plot the optimal proportional investment in the liquid asset 11(0, •), 
when p is respectively 0, -0.5 and 0.5. The Merton line corresponds to the case when I 
may be traded continuously, while in the case A = the proportion invested in I is actually 
lost. Notice that when A increases, 11^(0, z) goes from IIo(0, z) to 11^(0, z), increasingly or 
decreasingly depending on the value of p. 

Finally, it would be interesting to look at the impact of the parameters on the investor's 
behavior for t > 0. Analyzing that is more complicated: this cannot be performed by a 



35 




0.2 0.4 0.6 0.8 1 

A = - A = 5 A = 50 - 

A = 1 - A - 10 Merton 



Figure 4: Optimal investment (proportion) in the liquid asset II* (0, ■) as function of z with p = 0, 7 = 0. 

simple comparison of the feedback maps, since, as we have observed, these maps do not 
read the same variable as input for different values of the parameters. To overcome this 
problem one would need to perform Monte-Carlo simulations to study the distributions of 
the optimal investor's portfolio and strategies. We observe also that, since our (auxiliary) 
control problem is not autonomous, we cannot look at the stationary distribution as in 
[I]. However such approach would be intensive numerically, and for simplicity here we just 
consider the dependence of the feedback maps on the extra observation 16 ! for different 
values of 7, which still allows to illustrate the effect of partial observation. 

Since we are interested mainly in the impact of 7 on the strategies, we fix the other 
parameters, e.g. fix A = 5, p = 0. We consider an agent who at time t = 1 has a 
liquid wealth composed of X\ = 0.5 units, while having invested ao = Yq = 0.5 units in 
illiquid wealth at time (clearly, we assume that t\ has not occurred yet). We plot the 
optimal consumption rate and the optimal investment in the illiquid asset as function of 
the additional information, i.e. of B^\ which determines, together with 7, the value of 
Yi r*\ To be more explicit, we have to compute C*(l,Xi,Y\_), H*(l,Xi,Y\_) in terms of B^ 
and 7. Since from (|23p we have 

Yl = y oe ^- fl ^-hA=?7Bj 1) ) (91) 

16 We remark that in our model, what one really observes are the the processes L,I^\ so it would make 
sense to study the feedback maps as functions of them. However, our main interest is to study the dependence 
of the feedback strategy on 7. Since I^ 1 ' is still dependent on 7, to avoid the same problem pointed out 
above, we study their dependence on B^K 

17 This is true since p — 0, so that Y does not depend on W. 
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Figure 5: Optimal investment (proportion) in the liquid asset 11(0, •) as function of z with p = —.5, 7 = 0. 



by substitution we get the function to plot. 

In Figure [7] we observe the consumption rate: we see that the optimal consumption 
is increasing in which illustrates the unsurprising fact that when the agent knows 

that his illiquid investment is doing well, he should consume more. This effect is stronger, 
the more information the agent on the illiquid asset has (i.e., for 7 close to 1). Note that 
the impact of 7 here is higher than at time 0, as expected, since when time passes the 
importance of the additional information should increase. 

Figure [8] shows the optimal amount to invest in the risky asset. Again we observe the 
same behaviour with respect to -E^ 1 ) and 7. 
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2.5 



2 - 

1.5 - / / 

1 - / / 

0.5 - . 

' 1 1 1 1 1 

-3-2-10123 

7 = 7 = 0.4 7 = 0.8 

7 = 0.2 7 = 0.6 7 = 1 --- - 

Figure 8: Optimal proportion of liquid wealth to invest in the liquid asset in function of B[ for various 
7 (setting A = 5, p = 0, X% = .5, Yq = .5). 
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A Appendix 

Lemma A.l Given r > 0, for any (c, 7r, a) G -4.(V), there exists (c°,7r°) G .4o( r — a o) such 
that 

(c,vr)l {t < Tl} = (cV°)l {t < Tl} , dF®dsa.e.. (92) 

Proof. First, using the definition of G, by a simple monotone class argument, for any 
(c, 7r) which is (^)t>o-predictable we may find (c°,7T ) which is (W t V B t (1) ) t > -predictable 
satisfying ([92]) . It is straightforward to see that the admissibility constraint (c, ir, a) G .A(r) 
implies (c°,7r°) G >4o( r — ao), and the proof is complete. □ 

Proposition A. 2 

(i) G 7 is well defined on the set of measurable functions having at most linear growth. 

(ii) G 7 is linear and positive, in the sense that it maps positive functions into positive 
ones. As a consequence G 7 is increasing in the sense that 

<j> < $ G 7 [0] < G 7 [# 

(Hi) G 7 maps increasing functions in functions which are increasing with respect to both 
x and y. 

(iv) G 7 maps concave functions in functions which are concave with respect to (x,y). 

(v) Ifip(r) = r p , p G (0, 1), then (kj tP is defined in ([33]) ) 

G 7 [T/#,£x,£y) = ^G 7 [^](t,x,y), yt > 0, y(x,y) G M+, V£ > 0. (93) 
< G 7 [V#,x,y) < e kj ^{x + y) p , \/t > 0, V(x, y) G (94) 

(vi) Let p G (0, 1], and ip a p-Holder continuous function on R + . Then for all t > 0, 
x,x',y,y' > 0, and < h < 1, there exists some constant C > such that 

\G^mt,x,y)-G^[m,x',y)\ < C\x-x'\p, (95) 
\G^mt, X ,y)-G^mt,x,y')\ < Ce^y-yT, (96) 
\G-<li;}(t,x,y)-G~ f im + h,x,y)\ < C ie k ^yP h P/ 2 , (97) 

Proof. See [7] for the case 7 = 0. The general case is completely analogous, so we omit it 
for brevity. □ 

Lemma A. 3 Letp£ (0,1) and kL,Y,p> &J,p defined as in ([32]) - ([33]) . For (t, x,y) G [0, +oo), 

(C,7r) G A t {x), 



E 



^t,x,c,n + Y^y] < e^p^ix + yY, (98) 



for all s >t. In particular, combining (|98|) with Proposition \A .21 v) and denoting ip(r) = r v , 



E 



G 7 [( / 9](s,X*' a; ' c ' w ,y;*' s/ )l < e kj ^ kL ' Y ^ +kj ^ s -^(x + yf. (99) 
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Proof. See [7] for the case 7 = 0. The general case is completely analogous, so we omit it 
for brevity. □ 



Lemma A. 4 Set 



P( l -P)(„.2 2 , 2 2 



f(u L ,ui) := p(u L b L + ujbi) - (u L a L + UjGj + 2pu L u I a L a I ). 



Recalling (I12D . we have then k p = sup ltig ]R j , il/g [Q il ] f(ui, uj). For any b' Y , b'j such that 
b' Y + b'j = bj , define 

f b ' Y (uL,uy) := p(u L bL + uyb' Y ) - — V \ u 2 L cr 2 L + u Y aj(p 2 + 7 2 (1 - p 2 )) + 2puLUya L ai). 

f b ,(uj) := P6W-^; 2 (1-P 2 )(1-7V;, 

and k L,Y,P := su P«LeK,«y6[o,i] h' Y {uL,u Y ), k' Jp := f Vj {uj). Then k p < k' LYp + k' Jp , and this 
inequality is an equality if we choose 

b'y = 7 2 ^ + (l-7 2 )^. (100) 

Proof. Since fy Y (uL,uj) + fy } {ui) = f(uL,ui), then from the definition of k p ,k' LYp ,k' Jp 
we have 

k P = sup {fb'(uL,ui) + f b >(u I )) < k' LYp + k' Jp . 

« L GK,W/G[0,1] 

Since the maximizers of f-,fy Y ,fy J always exist, the inequality above becomes equality if 
and only if there exist a maximizer (u* L ,u Y ) of fy and a maximizer u*j of fy such that 
u Y = Uj. 

In the case 7 € (0, 1), by strict convexity of /// and fy , these maximizers are unique 
and may be computed explicitly with the first-order condition obtaining 

U} = PrOJ[ °' 1] ((l-pW(l- J p 2 )(l-7 2 ))' 



u* Y = Proj [01] 



b' 



( i_ p)(T 2 (1 _ p 2 )7 2y< 

Using that b' Y + b'j = &/, (jlOOp may be rewritten as 



b'j by 



(1 -7 2 ) 7 2 ' 

which implies u*j = u Y . 

To conclude, it remains to notice that under (jlOOp . for 7 = (resp. 7 = 1), fy does 
not depend on uy (resp. fy does not depend on uj), so that clearly in these cases we may 
choose u Y = u*j. □ 

Given (i, z) £ M + x (0, +00) and e G (0, z), we denote 

D £ (i,z) := [t,i+e) x (z-e,z + e) C R + x (0,+oo). (101) 
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Lemma A. 5 Let (t, z) G 1R + x (0, +00) and e G (0, z). 
1. There exist N E > such that for each (t,z) G D E (t,z 

$(t + h,z) -<&(t,z) 



lim sup 



h 



< N e . 



(102) 



2. •) G C l {(z — e, z + e);R) /or every t G + e) and there exist m E , M £ > such 
that 



m £ < $ z (t,z) < M e , V(t, z) G D £ (t, z). 



(103) 



3. $(i, •) is twice differentiate a.e. in {z — e,z + e) for every t G [t, i + e). Moreover, 
denoting by Of C (z — e,z + e) the set where $(t, •) is toice differentiable, there exists 
5 e > swc/i i/iai 



< -* e , Vt G [t,t + e), z G a 



Proof. 1. Setting 

J(t,x,y;c,ir) = E 



-(P+X)s 



U(c s ) + AG 7 [F](i + s, X° s ' x > c >", Y S °A ds 



since the equations for X, Y are autonomous, we have 

V(t,x,y) = sup J{t,x,y;c,7r). 



(104) 



(105) 



(106) 



Since we are in the power case, we have G 7 [y](i, x, y) = <3?° K(x + yJt) p - Applying Dynkin's 
formula to &°(x + yJt) p ', we see that G 7 [V](-, x, y) is differentiable and we get the estimate 



^[V](t,x,y) 



< C JiP (P[V]{t,x,y), 



(107) 



where 



1 



Cj, P = \bj\p + ^p(l-p)a 2 j. 



So we can differentiate (jlQ5|) with respect to t and, using (|107p and the nonnegativity of 
[/, we then get 



— J(t,x,y;c, tt) 



< C JiP E 

< C JiP F(t,x,y). 

Estimate (I108P is uniform in (c, 7r) G „4o(x), so from (|106p and the fact that 

|F(t + /i,x,y) - F(t,x,y)| < sup + h,x,y;c,ir) - J(t,x,y;c,ir)\, 

{c,ir)eAo(x) 



(108) 



we get the claim with 



N e = C JtP - sup V(t,x,y). 

D e (t,z) 
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2. Let (t, z) G D e (i, z). Since <£>(t, •) is concave, there exist the left and right derivatives 
z), $f(t,z), and <&~{t,z) > $+(t,z). To show that $(*,-) is differentiable at z we 
then must prove that the above inequality is actually an equality. Suppose by contradiction 
that 3>j(i, z) > $z(t, z). Consider the function defined for z\ G [z — e, z + e), ti G [t,i+ e) 

( p S (t 1 ,Z 1 ) = *(t, Z) + ^ Z) ± *» ( ' ! Z) (Z! - Z) - - Z) 2 + (JV e + ~ t). 

Due to item 1, the function <fr — <p s has a local maximum at (t, z) in (z — e, z + e) x [t, i+e) 
for each 5 > 0. Therefore the subsolution viscosity property at (t, z) implies 

-N.- t + *»•(«, ,) - % EM1£M _ A4 o r((> z) + ^ z2 1 

2 2 o 



rrl { $7(t,z) + $j(t,z) 1 _ ~\ 

sup F ct , ,--;c,0 < 0, 

■o.Sgm. V z / 



c>0, 06 

Letting 5 — > we get a contradiction as -y^ 2 ^ — > +oo, so ■) is differentiable at each 
z G (z — e,z + e) and for every t G [i,i + e). The fact that <&(t, •) G C 1 ((^ - e, z + e); R) 
for every t G [t, i + e) follows from concavity. Finally let us show (|103p . Let 5 = ^+^. By 
concavity of $(i, •) we have for every z G (z — e, z + e) and t G [t, i + e) 

*(t,z-e)-${t,z-e-6) < $ f < $(t, z + e + 5) - $(t, z + e) 

Since •) is strictly increasing for each t G [i, i + e) (Proposition 14. 1 1 ) . we have 

$(i,z- e) - ®(t,z - e- S) $(t,z + e + 5)-$(t,z + e) 

— ^ r - Lj < oo, — — — > 0. 



Calling 

&(t,z-e)-&(t,z-e-8) ., $(t,z + e + S) - $(t,z + e) 
M £ = sup , m £ = mf , 

te[i,t+e) d te[t,t+e) d 

by continuity of <3>, we have < m £ < M £ < oo, so the claim follows by (|109p . 

3. Let (t, z) G D £ (i, z). The fact that there exists a set Of with full Lebesgue measure 
such that <&(i, •) is differentiable at the points of Of follows from concavity of &(t, •) and 
Alexandrov's Theorem. Suppose that z G Of. Let 5,5\ > and consider the function 
defined for z\ G (z — e, z + e), t\ G [t, i + e) 

<p 5 (h, Zl ) = $(t, z) + $ z {t, z)(zi -z) + -(* w (t, z) - 8)(zi - zf - (N £ + dx)(h - t). 

Due to item 1, the function $ — ip s has a local minimum at (t, z) in (z — £,z + e) x [t, i+e) 
for each 5 > 0. Therefore the supersolution viscosity property at (t, z) and item 2 imply 

K 2 

N £ + 5 1 + # A $(i, z) - K 3 ^ £ - \<S>°r<(t, z) - -^-z 2 ($ zz (t, z) - 5) 

1 X 2 in 2 

-w2*P5riw 2 a (110> 
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Notice that, given ao, &o > and cq G R, there exists Qo > such that 

aoC"^ < c , £ < £ < -a - (111) 

Since we know that $ M < 0, from (jllOp we have that (lllip holds for £ = 2) — 5. So, 

we get the existence of 5 £ > independent of (t, z) G D £ (i, z) and of (5 such that 

*zz(t,z) < 5-S s . 

By arbitrariness of 5 we have the claim. □ 

Proposition A. 6 & is a viscosity solution in D E (t,z) of 

K 2 

-ip t + K x v-K 3 zp z -\f(t,z)<S> -^-z 2 (p zz - sup Hl v (<p z ,<p zz ;c,e) = 0, (112) 

c€[O,c],0e[-0,0] 



u>/iere 



eh K%5 £ 



Proof. The fact that $ is a supersolution of (|112p in D £ (i, z) is a straightforward conse- 
quence of the fact that it is a supersolution of ()66|) . as the supremum is taken over a small 
set in (|112p . Let us show that it is a subsolution in D £ (t,z). Take (t, z) G D £ (t,z) and 
let G C^A^i^R) be such that ip(t,z) = <S>(t,z) and in ip > $ in D £ (t,z). Since $ 
is once differentiable with respect to z, it must be (p z (t,z) = §> z (t,z). Now, if ip zz < —6 £ , 
then 

sup Hl v (ip z ,y zz ;c,6) = sup H^ v (ip z ,ip zz ;c,9), (113) 

c>o, eeR se[o,c], de[-e,0] 

so we have the desired subsolution inequality. Otherwise, suppose that ip zz (t, z) > —5 £ and 
consider the function defined for z\ G (z — e, z + e), t\ G [t, i + e) 

{p{ti,zi) = ip(t u z) + $ z {t 1 ,z)(z 1 -z)-^5 £ {z 1 -z) 2 . (114) 

We have 

<p(t u z) > (p(h,z) > ®(t u z), VtiG[i,t + e). (115) 
Now fix t\ G [t,t + e). Consider, for z\ G (z — e, z + e), the Dini derivative of 3> z at z\\ 

n+(h ,. v ,. ^(tl,Zl + /l) -^(tl,Zj) 

DJ$ z (ti,zi) := hmsup . 

Since $(£1, •) is concave, we have 

I>+$,(ti, *i) < 0, Vzi£(z-E,z + e). (116) 
Moreover, by Lemma IA.5I T3) . we have 

D+Q&uz!) < -5 e , Vz 1 eOt 1 . (117) 
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As consequence of (|1 16[) - (jl 1T() . of the fact that has full measure and of Lemma 3.3 in 
[9], we get by integrating two times (|117p 

< ®(t u z) + $ z (t 1 ,z)(z 1 -z)--S e (z 1 -z) 2 . (118) 
Combining ([TTSjl with (flT^ - pTSl) we get 

(p(t, z) = ®(t,z), (p>§ in (z - e, z + e) x [t,i + e). 
Now, since $ is a viscosity subsolution of (|66p . we have 

K 2 

— fit + K\ip — K^zifz — A < l )0 / 7 (t, z) ^-z 2 (p zz — sup Hi(<p z ,<p zz ;c,e) < 0, (119) 

On the other hand, we have 

sup Hl v (Cp z ,Cp zz ;c,6) = sup H^ v (ip z ,ip zz ;c,9), 

c>o,eeR £e[o,c], ee[~e,e] 

so that also 

if 2 

- (p t + if A <^ - tf 3 *0, - X^°F(t,z) - -fz 2 ^ zz - sup Hl(<p z ,<p zz ;c,0) < 0.(120) 

2 ae[o,c], ee[-e,e] 

Noticing that 

2) = <p(t,z), ¥^(M) = <p z (t,z), (p zz (t,z) > -S £ = <p zz (t,z), (121) 

and taking into account that H^ v is nondecreasing in the second derivative, combining (11201) 
and (|12ip we get the desired subsolution inequality for ip and the proof is complete. □ 

Lemma A. 7 Let a < b and F : [0, T) x (a, b) x E 3 — > R, (t,x,r,p,q) H> F(t,x,r,p,q) 
be continuous, continuously differentiable in (x,r,p,q), and proper degenerate elliptic (i.e. 
nondecreasing in r and nonincreasing in q). Let u € C 1,2 ([0, T) x (a,6);IR) be a classical 
solution in [0, T) x (a, b) to the parabolic equation 

u t + F(t,x,u,u x ,u xx ) = 0. (122) 

Then the space derivative v := u x is a viscosity solution in [0,T) x (a, b) to the parabolic 
equation 

v t + S/F(t,x,u(t,x),v,v x ) ■ (l,v,v x ,v xx ) = 0, (123) 

where VF = (F x , F r , F p , F q ) is the gradient of F(t,-). 

Proof. For x £ (a,b) and sufficiently small h > 0, define u h (t,x) = u(t,x + h) and 
v h = - . Then, due to continuous differentiability of u, we have that v h — > v locally 
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uniformly in [0,T) x (a, b) when h goes to 0. Furthermore, since u is a solution to (|122p . 
and using the differentiability of F, we see that 

v% = j^[ F (t,x,u h ,u x ,u xx ) - F(t,x,u,u x ,u xx fj 

= (vF(t M v h ,v h x )+£ h {t,x)) •(l,u h ,t£,t&), 

where 

£ h {t,x) 

:= / VF(t, x + hs, (1 - s)u h (t, x) + su(t, x),(l- s)v%{t, x) + su x (t, x), (1 - s)u h xx {t, x) + su xx (t 
Jo 

-F(t,x,u(t,x),v h (t,x),v^(t,x)). 

By continuity of F and of u, u x , u xx , and by the fact that v h (resp. v x ) goes to u x (resp. 
u xx ) as h goes to 0, we see that £ h goes to as h — > 0, locally uniformly in [0, T) x (a, b). 

So, it just remains to apply the stability result for viscosity solutions (see, e.g., [25] 
Prop. 5.9, Ch.4]), to get that v is a viscosity solution to the limiting equation (|123p . □ 
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